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Abstract 

In this paper we rigorously prove the existance of a non-trivial 
periodic orbit for the non-linear DDE: x'(t) = K s'm(x(t — 1)) for 
K = 1.6. We show that the equations on the Fourier equations have 
a solution by computing the local Brower degree. This degree can 
be computed by using a homotopy which validity can be checked by 
checking a finite number of inequalities. Checking these inequalities is 
done by a computer program. 



1 Introduction 

We will show that the equation: 

x'(t) = -Ksm{x{t-l)) (1) 

has a periodic solution for K = 1.6. By a solution we mean a C 1 function 
that satisfies the equation. 

Numerical simulations suggests that for | < K < 5.1 there is a attracting 
periodic orbit oscillating around 0. An informal argument for the existence of 
such an orbit also comes from the analysis of eigenvalues of the linearization 
[Wischert et al 94] . We will prove rigorously that for K — 1.6 this orbit 
exists. 

In the remainder of the introduction we briefly outline our method. We 
will use the Fourier coefficients of the periodic orbit. First we rescale the 
time to obtain a 2ir period. We make a substitution: x{t) = x (-) where r is 
a parameter. If r is equal to Ipr (where T is the period of the solution in the 
original equation) we will have a 27r-periodic solution of the new equation: 
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-// N 1 / ft\ K . ( ft 
x (t) = —x I — = sin I x 



T \T J T \ \T 

K . ( (t-r\\ K . 
= sin I x I I I = sm(x(t — r)j 

By renaming x to x we obtain: 

K 

x\t) = sin(x(t-r)) (2) 

T 

We will prove for this equation that there exist a r from a small interval 
such there exist a 27r-periodic orbit in a small neighborhood of a specified 
function. Note that we don't obtain the exact value of the period (even if 
numerical solution suggests it is 4, i.e. r = |) but treat r as a variable. 

To prove the existence of the orbit we will use the method of self-consistent 
bounds that was introduced in the context of Kuramoto-Shivashinsky PDEs 
in |Zgliczyhski 01| , |Zgliczyhski 04| . When applied to the boundary value 



problem for ODEs or DDEs this method is similar to the Cesari method 
introduced in [Cesari 64] but doesn't require one of the conditions - see 
Section 2.4 in |Zgliczyhski 01| for a comparison. Below we briefly describe 
the method. 

We will derive from (j2J) the equation on Fourier coefficients. Let's first 
write ((2j) with a Taylor expansion instead of sinus: 



k 



ra=0 v ' 



2k+l 



The equations on the Fourier coefficients will look similar 

t f-^ (2fc+l)! v Jn 



inr 



k=0 

where (c*^ 2fc+1 ^)„ means convoluting the sequence c with itself 2k + 1 times 
and then taking the n-th coefficient. The operation of convolution and this 
notation is introduced in details in section El 

We will show that such a function has a non-trivial zero: 



F{t,c) = |mre inr c, 
Where the domain is: 



+ kY" [ ' (c* (2k+1) ) I 
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X, 



Vn G Z : \x„ I < 







(H + i 



for some /3 > 0. On we consider the product topology (which is equivalent 
to component-wise convergence). It turns out that on such a domain the sums 
in the definition of F are convergent, F is continuous and each sequence c G 
X/3 corresponds to a real-valued continuous function. Moreover if F(c) = 
then the function is C 1 . 

We will search for a zero in the neighborhood of an approximated solution 
obtained from numerical simulations. We will denote by c this approxima- 
tion. For n < n < 5 it is equal to: 



n 





-0.1521000000 - 0.1163508047i 


0.0001123121 - 0.0002746107* 


-0.0000008173 - 0.0000001014? 



For n > 5 the c n is zero, for n < we have c n = c_ n . The r in the 
approximation is: f = 1.570796. 

For each I G N let us define the Galerkin projection Pi and immersion Qf. 



P l :£ L 3c^ (co, . . . , Cj) 6 M X C 1 

Q r .RxC l 3 (c ,...,q) -> (...,0,0,c=r,...,co,...,Q,0,0,...) ec z 

Let's note that we need only the non-negative terms in the finite space 
as the negative terms can be obtained by conjugation. Also as cq = c~o so 
we have Co G R. From the compactness of Xp it will be easy to show that: 

Lemma. \15\ Let (3 > 0, lo > 0, r, r G E be fixed. If for each I > lo there 
is a c l G X/3 and Ti G [r, r] smc/j t/iat PiF(ti,c 1 ) = iften t/tere eMsi a 
(c°,r ) G X/3 x [r,r] swcft ttai F(r ,c°) = 0. 

Moreover if all the c 1 are in a closed set D then c° G D . 

Thus it's enough that there is a zero for every Galerkin projection of F: 

Fi : K x (l x C 1 ) 9 (r, c) -> PiF(r, Q,(c)) GlxC 1 

This will allow us to use a finite-dimensional topological method to show 
a zero of F. Let's note that the real dimension of the domain is 21 + 2 
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while of the image is 21 + 1 so if there is a zero then one can expect a 1- 
dimensional manifold of zeros. This manifold can be easily identified - if x(-) 
is a solution then x(- + 0) is also. This means that if Fi(t, Co, ci, . . . , q) = 
then Fi(t, cq, e^ci, . . . , e tl ^ci) = what can be easy checked. To use the 
topological method we want to have the zero isolated so we will limit the 
domain - we assume that c\ G C\ + R. Of course finding a zero in a limited 
domain implies a zero of the full system. By F t we will denote the Fi with 
the smaller domain. 

The before-mentioned method to prove the existence of a zero of Fi is 
the local Brouwer degree (introduced e.g. in [Smoller 83] ). Let's denote by 
deg(-F}, U,x) the degree of x on U. It is known that if the degree is non-zero 
then 3y G U : f(y) — x. We will use a neighborhood of c as U and x — 0. 

To compute the degree we will use the homotopy invariance of the local 
Brouwer degree. As our homotopy Hi we will use a linear deformation of Fi 
into a function Gi that contains the most important terms of Fi (that is not 
strictly a linearization of Fi but it is close to): 

H(h, r, c) = hF(r, c) + (1 - h)Gi(T, c) 

The degree of G\ will be easy to show that it is non-zero. What we will 
need to show is that ^ Hi([0] 1]; dU). To show that it's enough to show that 
the terms in Gi dominates the (mainly non-linear) terms that are in Fi but 
not Gi - that \Gi\ > \Fi — Gi\. That part of the proof is computer-assisted - 
the proofs of the estimates are in this paper but computing the exact values 
and checking that the inequalities holds is done by a computer program using 
the CAPD package for rigorous interval arithmetics. 



1.1 The program 

This program is written in C++ and can be downloaded from 



http://www.im.uj.edu.pl/MikolajZalewski/dl/delay-sin.tgz . The rounding- 



mode changing code required by the interval arithmetic is system-dependent 
and has been checked to work on PCs (both 32-bit and 64-bit) on both 
Windows (compiled with cygwin) and Linux (compiled with gcc). It should 
also work on SPARC and Mac OS X although that has not been tested. 
Using other CPUs or compilers might require modifications to the rounding 
code. 
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2 Fourier coefficients 



We will use the following notation: if {x n }'^L_ 00 and {yn}^L-oo are sequences 
with complex values then by x * y we will denote the convolution of the 
sequences: 



x*y:= I ^ x k y n - k 



-oo 



It is well known that if the sum converges the operation of convolution is 
associative. We will also use the notation (...)„ for the n-th coefficient of the 
sequence in brackets — e.g. {x * y) n = ^^=—00 x kUn—k- 

Also: 



x* := p * ■ ■ ■ * x, (for k > 1) 
k times 

1 if n = 



( x * 1 Qifn^O 

Of course the sum in the definition of the convolution may be not con- 
verging so we will limit our attention to a domain where the sum will be 
always convergent. The following lemma holds: 

Lemma 1. If x and y is such that 3a > 2,Pi,{3 2 > : Vn : \x n \ < 
(| ra |+i) a ' \Vn\ - (R+IT ' tflen is convergent for each n and \ (x*y) n \ < 

Cj JM^ where C = 222+11 

(M+l) a— 1 

Proof: Let's assume n > and let's estimate \{x * y) n \: 



00 




00 




n—1 




00 


^ ^ ^kVn—k 


< 


X-kUn+k 


+ 


^ ^ %kyn—k 


+ 


%n+ky~k 


k=— 00 




k=0 




k=l 




k=0 



En— 1 <~ R R 1 1 1 

k=l X kUn-k\ ^ PlP2 2^=1 (fc+i)<* ( n _fe+i)a 

— 2 A/^2 X!l<fc<§ (fc+l)« (n-fc+l) a 

— 2 A/^2 Za<fc<f (fc+l)" ((n/2)+l)<* 
^ o /3i/3 2 2" V- 1 

— Z (n+2)« Z^l<fe<f 
^ p /3ife2" V- 1 

— ^(n+2)« Z^fc>l (fc+l)« 

Efc=0 X -fc^n+fel — A/^2 Sfc=0 (fc+l) a (n+fc) a +l 
— ( n +l)« Z^fe>0 (fc+l) Q 
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Analogically: 



fc=0 



< 



r— 



{n+l) a ^{k + l) a 



The sums 



OO 

k=k k a +l 
oo 



q-j- can be estimated by integrals: 

1 



y^— <i 



fc=0 



{x + l) a 



dx 



a 



00 1 1 ^00 1 



a — 1 
a + 1 



^—dx - 
l) a 2 a (a-l) 



From that we obtain: 



< 2(2a + l) 

~~ (n + l) a a - 1 

The result for n < can be obtained by analogous estimations or by 
taking a sequences x, y: x n := x_ n , y n : = y_ n and applying for them the 
result for n > 0. 

□ 

With one exception we will use this lemma for a = 2. For a = 2 we have 
C = 10. 

Observation 1. If x is such that 3a > 2, (3 > : Vn : |x n | < Q^mys 
I (***)»! < ^(RTTT w/iere C = 2J2 ^T 

Lemma 2. // x and y are such that x_„ = x^, y_ n = y^ then (x * y)- n = 
* V)n 

Thus we have that if x G Xp 1 and y G -X/? 2 then x * y & Xcp^ (where 
X/3 was defined in the introduction and C is from Lemma CD). It will be also 
useful to define a set of sequences from any X@: 

X:=\JX P 

We will limit ourselves to the sequences in X. For them we have: 

Lemma 3. Let c G X. Then X]^L-oo c n^ mt converges to a real-valued con- 
tinuous function 



< 



(71 + 1 



a + 1 
2<*(a - 1) 
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Proof: The functions ^2 n= _ N c n e mt are continuous and real- valued 
as c n = ~c~^ =>■ c n e mt + c_ n e~ m * G R. They uniformly convergent as 

Zl^L-oo l c n em *l < Z]^L-oo (|n|+i) 2 < 00 - Thus tne l imit is a ^ s0 real-valued 
and continuous. □ 
Let's note that c G X doesn't guarantee that the function is C 1 : 
We use the convolution because of this well known fact: 

Lemma 4. If c G X are the Fourier coefficients of x(t), d G X are the 
coefficients of y(t) then the Fourier coefficients of x(t) ■ y(t) are c* d. As a 
consequence the coefficients of x n (t) are c* n . 

Now we can write the equation j2]) on the Fourier coefficients. 

Theorem 1. Let r be fixed and x{t) : R — > R be a 2n -periodic function with 
the Fourier coefficients c G X . Then: 

(i) for each n G Z the sum on the right-hand side of the following equation 
converges: 

k=0 K ' 

(ii) x(t) is a 2ir -periodic solution of (Q|) 4=> c satisfies the equations 

To prove the theorem we will need some lemmas. As it has been already 
mentioned that c G X doesn't imply that x{t) is C 1 . However if c is the 
solution of equation Q then we have the following lemma that will allow us 
to show that x(t) is C 1 (x(t) is even C°°): 

Lemma 5. If c satisfies equation (GjJ and (3 > 0, a > 2 are such that Wn : 

\ Cn \ — (H+TF ^ ' ^ U ' \ Cn ^ — (M+1) Q+1 

Proof: Let C := *2s±*l. We have that: 

a— 1 



A 



(-i) fc 

k=0 (2fc+l)! 



(c*(*+D) , 



We also have \inc n \ 
equal hence we obtain: 

K 



n c„ 



< 
< 



K 



Eoo 1 
fc=0 (2fc+l)! 



( c *(2fc+l)^ 

2k I3 2k+1 



C 



K 



k=0 (2fc+l)!^ (\n\+iy 

1 - sinh(C/3) 



(4) 



r C(|n| + l) a 



t C(\n\+1)° 

The two sides of the equation must be 
sinh(C/3) > \n\ \c n \ 
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Thus for n 7^ the /3' = ^ sinh(C/5) satisfies the assertion. If it is not 
satisfied for n = we can increase the (3' . □ 
Increasing the a is important as we have: 

Lemma 6. Let c be a sequence of complex values satisfying |c„| < ^ n \ + \yi 
for some (3 and c n = c_ n . Then the sequence c is a sequence of Fourier 
coefficients of a real-valued C 1 function x{t). 

Proof: The sequences Ylk=-n c k^ lkt and Ylk=-n ^ c fe e * fc * are real- valued, 
the second is the derivative of the first one and are uniformly convergent as 
n — > oo. Hence both converge to continuous function and the second is the 
derivative of the first one. Thus x(t) = YlT=-oo c k^ kt is C 1 . □ 

Proof of the Theorem Q} 

Ad (i): The convolutions are convergent because x G X. From the 
equation ([1]) from the proof of Lemma [5] we have that if x G Xp then 



Ad (ii): Implication =^: It's enough to show that {mc n }^l_ 00 are the 

(2fc+l)! 



Fourier coefficients of x'(t) while -f E^olfeS)! (c< 2k+ ^) n e~ inT are the 



Fourier coefficients of — — sina;(t — r). 

The first part can be obtained by integrating J 27r x'{t)e lkt dt by parts. 

As for the second we have that — ^ Y^k=o (2fc+i)! — r)) 2fc+1 — > — y sinx(t- 
r) as JV — > oo. The Fourier coefficients of x(t — r) are equal to d :— 
{c n e~ mT } c ^ = _ 00 . Hence the coefficients of x(t — r) 2fc+1 are cH 2fc+1 ) which is 
equal to c *( 2fe+1 ) e~ inT . 

Hence we have that the Fourier coefficients of — ^ ^^=0 (2k+iy. ( x (^ ~ T )) 2k+1 

are equal to j-f J2k=o (fe^yr ( c * {2k+1) ) n e ~ inT }°° ■ lt has been shown 
that this sequence in convergent as iV — ► 00 thus the n-th coefficient of 
-f smx(t - t) is -f J2T=o (ferSl (c*( 2fc+1 )) n e- inr . This ends the proof of 
this case. 

Implication •<=: from Lemmas [5] and [6] we obtain that x(t) is a C 1 function. 
We know that -f £~ ^=^y (c* (2fc+1) ) n e" inr are the Fourier coefficients of 

— — sina;(i — r). They are equal to inc n - the Fourier coefficients of x'(t). 
So both — — sina;(t — r) and x'(t) are continuous and 27r-periodic functions 
with equal Fourier coefficients. Hence the functions themselves are equal and 
equation j2|) is satisfied. □ 
Let's define a function F : R x X — * C z (with the product topology on 
C z ) corresponding to the equation (jHJ): 
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F(r,c) := Lr^c. + K±±^L(^)\ 

Of course having F(r, c) = is equivalent to the fact that r, c satisfies 
equation (J3j) . In the rest of the paper we will show that F has a nontrivial 
zero. Let's note that F(r,c) n = F(r, c)_ n . 

As we will use topological tools we will want F to be continuous. First 
let's note that the convolution is not continuous on the whole X (as written 
above we use the product topology on X) but we have: 

Lemma 7. The operation of convolution is continuous on each Xp 

Proof: Let's fix some x°,y° G Xp. Let's fix some 5 and take some 
x, y G Xp from a neighborhood of (x°,y°): x,y are such that V|n| < N$ : 
\x n — (x°)n\, \y n ~ (z/°)n| < S where Ns is large enough that for \n\ > N$ : 
x,y,x°,y° e Xp ^ \x n - (x°) n \, \y n - (y°) n \ < b~- We have that: 

\(x° *y° - x *y) n \ <\ (x° * (y° - y))J + | - x) * y) n \ 



oo co 1 

-° * (v° - »)).! < E | K),-! • | (A,-, - M £ & E (]7]TT 



J=-CO 



j=-oo 



This tends to zero as 5 — > 0. After an analogous estimation for | — jc) * y) n 
we have that the convolution is continuous. □ 



Lemma 8. The function F is continuous on each Xp. 

Proof: The convolutions are continuous on each Xp and the result of a 
convolution lays in a Xp for some (3' so K^2^ =0 ( c * ( ' 2fc+1 ' ) ) n i s contin- 

uous for each N. From estimations as in equation (j3J) we have that for each 
coefficient this series converges uniformly as N — > oo hence we obtain that 



the limit K 



(-l) fc (>(2fc+l) 



(c*( 2k+r> ) is continuous. 



/fe=0 (2fc+l)l 

The term inre mT c n is also continuous so F is continuous. 



□ 



3 Some estimates 

As mentioned in the introduction we will need some estimates to show that 
the inequality holds in the neighborhood of c. We will use two kinds of sets: 
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c + X/3 2 that will be used to obtain finer estimates and Xp 1 (where (3\ will 
be large enough to contain the whole set c + Xp 2 ) for some more rough but 
simpler ones. 

In this section we will only assume of c that almost all coefficients are 
equal to zero. By Y\ we will denote the space of the possible values of c - the 
set of sequences such that at most the elements — /, . . . , I are non-zero: 

Yi := {c G X : Mn : \n\ > I c n = 0} 
First let's note two simple properties: 

Lemma 9. If c G Y x then c* k G Y kl 

Lemma 10. Ifx,yEX then 

k=n 
k=0 

To estimate ^ ( x * p ) n ~ an element of the sum in the equation j3]) - for 
x G X/s it's enough to use LemmaHl However for the sets of the form c + Xp 
we will use a more sophisticated estimation: 

Lemma 11. Let c G Yi, p > 1, (5 > 0. Then for any x G Xp: 

P \ \ ^ I *k\ (-y P -k-i fl P \ 

Where C = 10 

Proof: From Lemmas [T1 and [TUI we have: 



K 

-ri 
pi 



c + xyj> 



< 



K 
p\ 



fc=0 ^ 




< K ( \r*P\ 4- Y^P -1 ( P | V^ 00 |_*fc| rip-k-l P p k \ 
- p\ \ \ C n I I" Z^fc=0 U I Z^j=-oo IS I U (\n-j\+l) 2 J 

We have c 6 so c* k = for j > kl what ends the proof. □ 
Of course we can use the previous lemma only for a finite number of terms. 
To estimate the tail we will use a weaker estimation by using a neighborhood 
of the second type, applying Lemma [Hand summing the geometric sequence: 



10 



Lemma 12. Let N be odd, (3 £ [0; £jA . If x £ i/zen 



A' 



, (2* + l)! 



x 



*(2fc+l) 



WTzere C = 10. 

Proof: Using Lemma [U 



< 



A ■ (CP) 



N-l 



Eoo 



K 



(2fc+l)! 



x *(2fc+l) 



< E 



< 



K 



k={N~l)/2 N\N 2k ~( N 
K s^oo (Cf3) 2k 



V 

AT! Zuk=(N-l)/2 N 2k ~( N - 1 ) (|n|+l) 2 



irl(** (a * +1) ), 

/3 



iV-l 



TV! 
K(C(3) N - 



(H+iy 



iV!(l_(^) 2 ) (|n|+l) 2 

The geometric sequence is convergent because /3<^=^(^) 2 <1. □ 
From Lemma [TT] we can obtain an estimation that after multiplication by 

(|n| + l) 2 is independent of n. That will be used to estimate the term in all 

inequalities for high n-s with one formula: 

Lemma 13. Let > 0, c £ Yi, p > 1, N > pi. Then for each x £ Xp, 
n > N we have: 




Proof: Note that n > pi 



and 



(|n-i|+l) a 



(| n _ i |+l)2 ( n+1 )2 



(Ar-j+ip ■ (^ttf as N > J and apply Lemma [n] □ 

For the first terms we will need to have a better estimation than in Lemma 
[TT1 so we will regroup the terms (to understand why this regrouping helps 
let's compare two estimates: |1.1 • x — 1 • x\ < l.l\x\ + l\x\ = 2.1\x\ and 
1 1.1 -x-l-x\ = |(1.1 - l)x\ = 0.1\x\). 

Lemma 14. Let x £ Xa, c £ Y it 



mre 



(Cn ~\~ X T: 



(-l) fc 

fc=0 (2fc+l)! 



V C + X 



*(2fc+l) 



(mre 



+ A) c n - 



§(c*\ + %(c*%-%(c*X + 
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Where lpJ = -f ( * ) + fi ( jj ) (<" 5 -»), - fr ( J ) (^""J, 



and we assume 




= /or k > n. 



4 Proving the existence of the periodic solution 

We will search for the orbit in the neighborhood of c and f defined in the 
introduction. Let's define the sets and boundaries on which we will operate: 

X\ := Xp 1 

X 2 := c + Xp 2 
X 3 := {yeX 2 :y 1 -c 1 E R} 
t '.— t — Ar 
r := f + Ar 

where Ar = 0.000001, & = 0.0000002438, A = 0.766763. 

We will prove that a solution exists in the set [r;r] x X 3 . The f3i, (5 2 are 
such that X 3 <Z X 2 <Z X 1 . 

As mentioned in the introduction we will use the Galerkin projections 
of F with the condition y-i — C\ 6l added to make the zero isolated. The 
Pi and Qi were defined in the introduction as the projection and immersion 
of the finite-dimensional space. Let's first prove the lemma stated in the 
introduction. 

Lemma 15. Let {3 > 0, l > 0, r, r G ffi be fixed. If for each I > l 
there is a c l G Xp and t\ G [r, r] smc/j £/ia£ PiF(ri,c l ) = £/ien £/iere existf 
(c°,r ) G x [r,r] such that F(t ,c°) = 0. 

Moreover, if all the c l are in a closed set D then c° e D. 

Proof: From the compactness of [r. r] x Xp there exists a subsequence Ik 
such that c lk converges to a limit c° and T\ k converges to To. Let's fix n > and 
note that for I > n if Pi(F(t, c)) = then F(t, c) n = 0, so from the continuity 
of F we have F(r,c°) n = lim*^ F(r lk , c lk ) n = lim ^oo Pi k F(n k , c lk ) n = 0. 
For n < we have F(t , c°) n = lim^^ F(n k , c lk ) n = lim^^ F(r Zfc , c^)_ n = 
Thus F(r , c°) = 0. 

The last assertion is obvious. □ 
Before defining the homotopy let's introduce some auxiliary notations. 
Let's denote: f n (j) '■= inre mr . By L n we will denote the linear part of f n : 
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L n (r) := (ine mT ° — n 2 e mT °){r — r ). By r n we will denote the non-linear part: 
r n (r) ■= f n (r) - / n (r ) - Ljr). 

We will define the homotopy on the whole infinite-dimensional space: 
H : [0; 1] x [r;r] x X 3 — > C z . We will use only the Galerkin projections 
of this homotopy to deform Fi (that is the Galerkin projection of F - see 
introduction) but to prove that there are no zeros on the boundaries it's a 
bit more convenient to define the homotopy on the whole space. As written 
in the introduction the homotopy is a linear deformation of F to a nearly 
linear function G: 

H(h, r, x) := hF(r, c + x) + (l- h)G(r, x) 
where G on the n-th coefficient is equal to: 



(f n (r)+K) x n forn^il 

(/nOo) + K + 71,0 + 7i,2n) %n + L{r)c n for n = ±1 



G(r,x) n := 

Where j Pi j is from Lemma | 

The G is not strictly a linearization of F as it doesn't contain all the 
linear terms and is not linear with respect to r. However it contains the 
most important terms - the rest will be shown to be small compared to them. 
Let's note that for n = ±1 we have x n _ 2 „ = = x n (because x±\ G M) so 
the term with 7i j 2n in F will be a linear term with respect to x n . 

One can write explicit formulas for H. For n ^ ±1 we have: 

H(h, r, x) n := (/ n (r) + K) x n + 
+h ((/ n (r) + K)c n + Er=i (^T)I ((* + c)* (2fc+1) ) n ) 

And for n = ±1: 

H(t,r,x) n := (/„(tq) + K + 7i )0 + 7i,2n)^n + -^„(r)c n + 

+^ ((/„(t ) + c n - f (c* 3 ) n + f (c* 5 ) n - f (c* 7 )J + 

+hr n (r)c n + h(f n (r) - f n (r ))x n + hR(c + x) 
Where R is a short notation for: 



oc 



R(y) ■= V - V* (2k+1) 

To show that any Galerkin projection of H doesn't have a zero on the 
boundary of [r;r] x Pi(X 3 ) it's enough to show that H doesn't have a zero 
on the boundary. More exactly: 
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ft 



(|n|+l)2 



Theorem 2. Ze£ x G X 3 , r G [t, r]. Let x be such that 3n : \x n 
or let re {t,t}. ThenVh G [0; 1] : H(h,T,x) ^ 0. 

Proof: The proof is computer assisted - the calculations are done by the 
program. The computation is as follow. 

First let's note that if for n < we have \x n \ = then this is also 

true for —n. Thus it's enough to consider n > (and r). We do it in two 
steps: 

1. Let's assume that \x n \ — (uj+ip for some n > 0, n ^ 1. We will prove 
that: 



K/nGr) +#) 



> 



{Ur) + K)c n + Y.t 



K 



((x + c)*( 2k+ V) n 



=1 (2fc+l)! 

As the value of H (h, r, x) n is the difference of these two terms (with the 
second one multiplied by a h < 1), it will follow that H(h,r,x) n ^ thus 
H(h,r,x) ^ 0. The LHS will be estimated by: 



\{f n (r) + K)x n \ > \\f n {r)\-K\ 



02 



(5) 



(|n| + l) 2 

and |/ti(t)| = ti\t\ is an interval in the interval arithmetic. 

The terms of RHS are estimated depending on n. For n < 225 we can 
estimate each n separately. For n > 225 we will want to check all the 
inequalities in some finite computations. 

Let's note that for n > 255 we have (f n ( T ) + K)c n = as c n = 0. In the 
sum YlkLi (2k+i)\ (( x + c)*^ 2fe+1 )) re we will estimate the terms k = 1, . . . , 22 
from Lemma [T3] while the term k > 22 from Lemma [T21 Now if we multiply 
the estimates for both LHS and RHS by (\n\ + l) 2 they will be independent 
of n. Hence we will check all the n > 225 by checking only one inequality. 

For n < 225 we have to use some more sophisticated estimations. Let's 
use Lemma [H] to group the terms: 

Un(T) + K)d n + £r=l J^y {(x + c)*^) 

((/n(r) + K) c n - f (c*% + - £ (c 

J2 P =i J2j=-ei 1p,j x n-j + R(c 
< \(f n (r)+K)c n -§(c*\ 

\(fn(r) - /n(r ))c n | + 
+ |(/ n (T )+K)c n -f (c* 3 )„ + f (c* 5 

v-7 V 6i I I fe(10fe) p - 1 
"i" Z^p=l Z^=-6Z I IPJ\ (\n-j |+1) 2 

+Ef =4 |((^+£r (2fc+1) )J + |E 



,*7 



)„) + 



< 



+ |i?(c + 



n 



+ 



;r 



+ 



-\oo 
;k=23 



((x + c)*( 2fc+1 )) ) 
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Table 1: Estimates of LHS and RHS for small n 



n 


LHS 


RHS 





4.1792 


10 


-7 


0.1687 


10 


_7 


o 

I 


0.4474 


10" 


-7 


0.0958 


io- 


-7 


3 


0.5080 


10" 


-7 


0.0827 


10" 


-7 


4 


0.4893 


io- 


-7 


0.0160 


io- 


-7 


5 


0.4537 


10" 


-7 


0.0118 


10" 


-7 


6 


0.4171 


10" 


-7 


0.0070 


10" 


-7 


7 


0.3834 


10" 


-7 


0.3465 


10" 


-7 


8 


0.3536 


10" 


-7 


0.0040 


10" 


-7 


9 


0.3274 


10" 


-7 


0.0036 


10" 


-7 



The estimates obtained by the program. The LHS is closest to the RHS for 
n = 7 because we have artificially chosen c-j — 0. That makes the term 
\{f n (T)+K)cn-§ (c* 3 )„ + f (c* 5 )„~ fr (c* 7 )J f or n = 7 equal to approximately 
0.3411 ■ 1CT 7 . 

in the first term of RHS we will use the estimation: 

\(fn(r) - fn(r ))\ < max re[li7] \f'(r)(r - r )| < 

< max re [ z>T ] \ine mT — n 2 re mT \ ■ |Ar| < (6) 

< max re [ r>T ] (n + n 2 r) | Ar | 

The fact that it grows quickly with n is irrelevant as for n > 5 we have 
c„ = 0. 

The | (/ n (r ) + AT) c n - f (c* 3 ) n + f (c* 5 ) n - f (c* 7 )„| can be directly com- 
puted and is small as it is the numerical solution that is close to zero. 

The e;=i e;=. -6* ITpj'I ^nn-ji+i) 2 can ^ e directly computed and is small 
for p > 1 as 10/?2 is small. If the terms for p = 1 wouldn't be small we could 
try to diagonalize linear part. But there is no such need - we check that 
the inequality holds. We estimate | ((x + c)*( 2fc+1) )J for k G {4, ..,22} from 
Lemma fTTl To estimate \ J2T=23 (i x + c)*^ 2fe+1 ^)J we use Lemma[l2l In both 
estimates we have ( 2 fc+i)! w ith k large enough to make the result small. 

Thus we see that unless the 7^. are big, the estimation of RHS should be 
small. The program computes the estimation for LHS and RHS, compares 
them and for each n in 0, 2, 3, 4, .., 224 the inequality to hold. This ends the 
case for n 7^ 1. 

2. We have two cases left: r G {z,t} and |xi| = We have two 
variables left but only one equation - the equation for n — 1. However this 
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is a complex equation and the variables are real so we will be able to show 
that in both cases H(r, x) ^ 0. 
Let's denote: 



L T := Li(t)ci 
L x ■= (/»(fo) + K + 71,0 + 71,2) JCi 

iV := (/ B (T ) (^ 3 )n+f V (^)n-f (£* 7 )n + 

+i2(c) + r n (r)c n + (/(r) - /(r ))x„ 

The L T ,L X ,N depend on r, x but to make the notation short we won't 
write it. They are chosen such that G(r,x)i = L T + L x and F(r,x)x = 
L T + L x + N and H(h, r, x)\ = L T + L x + hN. If r G {r, r} then we show 
that: 



|L r | > + \N\ 

The |L T | is estimated as: 



\L\(t)ci\ = \in — n \ ■ At ■ \c 



V2A 



T ■ c i 



71=1 



We estimate the terms in \N\ as for n 7^ 1 with the exception for the 
|(/(t) — /(r ))x n | and \r n {r)c n \ that were not present in the previous case. 
The the first one we use the inequality ©. For the second we use the 
estimation: 

\r n (r)c n \ < sup TG[z;T] \\f"{r){T - r ) 2 | ■ \c n \ 
= ^ sup^.-j |- (2n 2 + in 3 r) e inr \ 
= ^sup re[z .-] |2 + ir| 

Having all these estimates our program checks that the inequality is sat- 
isfied. 

In the case = 4 the inequality \L X \ > \L T \ + \N\ is of course false. To 
prove that there is no zero will need to use the fact that r Gl and X\ G JR.. 
On the figure 1 we have a sketch how the sets L T and L x + N looks like - we 
see that they shouldn't intersect. 

Formally we will show that {tan arg(L r ) : (t,x) G [t,t]xX 3 ,Xi = ±^}fl 
{tanarg(L x + N) : (r,x) G [r, r] x X 3 ,xi = = (where arg is the 

complex number argument). Of course this implies that the sum of these 
terms is different from zero. 
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Figure 1: The possible values of L T and L x + N 



The arg(L T ) is easy to compute as this is tanarg ((ie lT ° — e tT °) c\) and 
tanargz for a complex number z can be computed as ||. To estimate 
the other set we will use the estimation for N from the previous point - 
let's assume that \N\ < A for some A > 0. We have that L x + N G 
± U'n( T ) + K + 7i,o + 71,2) if + [ — A; A] + [— A; A]i and using the interval arith- 
metics we can find tanarg(L x + iV). The program checks that these two sets 
are disjoint and this ends the proof of the theorem. □ 

To finish the proof of the existence of the orbit let's define a second 
homotopy. 



H L (h,r,x) := 
It deforms G into: 



Uni.ro + h(r - r )) + K) x n for n^il 
(fn(r ) + K + 71,0 + 7i,2n) + £( r ) c i for n = ±1 



G i (/i,r,x) := 



(fn(ro) + K) x n for n^±l 

(/n(To) + ^ + 7l,0 + 7l,2n) x n + L(t)ci for n = ±1 



Lemma 16. If (r, x) G [r, t] x X 3 swc/i i/zai r G {r, r} or 3n : |x n | = (i^j+ip 
^en H L (h,x, r) ^ 
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Proof: Let h, x, r be such that H L (h, x, r) = 0. 

Let n 7^ ±1. The r + h{r — r ) G [r; r] thus |/„(r + /i(r — r )) + K| can 
be estimated as in equation ([5]). For each such value we have proven that it 
is strictly greater than an RHS> 0. Hence \f n ijo + h(r — r )) + K\ > and 
H L (h,r,x) = => x n = 0. 

But if x n = for each ?i / ±1 then H L (h,r,x) = G(t,x) = H(0,t,x) 
and has no zeros on the boundary. □ 

Observation 2. The Galerkin projection of H: i.e. Hi(r,yo, . . . ,yi) := 
PiH(r,Qi(y , . . .,y t )) (for (y , ...,y n ) G M x M x C n ~ l ) is a homotopy from 
F 1 (t,c+ ■) to the projection of G. Analogically the projection of H L is a 
homotopy from G to G L . There are no zeros on the boundaries for these 
homotopies. 

Observation 3. The degree deg(Fi, [r,r] x P/(X 3 ),0) is well defined. 

Proof: The homotopy Hi for h = 1 have no zeros on the boundary and 
this is F h □ 

Lemma 17. deg(F ; , [r,r] x P t (X 3 ),0) + 

Proof: We know that the degree of Fi is equal to the degree of the 
Galerkin projection of G L - let's denote it by Gf. The Gf is a linear function. 
If the determinant of the differential was zero then there would be a zero on 
the boundary of any neighborhood of (tq, 0) so the determinant is non-zero. 
That means that the degree is ±1 thus it is non-zero. □ 

Theorem 3. The equation (dp has a periodic solution for some r G [r, r] 
whose Fourier coefficients are in the set X 3 . 

Proof: From the Lemma [17| we know that the local Brower degree is 
non-zero thus each Fi has a zero in X 3 . From Lemma [15] we have that F has 
a zero in X 3 . From Theorem Q] we obtain that in X 3 there is a solution of 
the equation. □ 



5 Conclusions 

In this paper I was able to rigorously prove the existance of a periodic orbit 
for K = 1.6. I was not able to show that the period is 4 even if the numerical 
simulations suggests that - in the proof we need the r as a variable for the 
image and the domain to have the same dimension. 

The value K — 1.6 has been chosen because it's easiest to find the c n 
values for an attracting orbit. However using the Newton method it should 
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be possible to find an approximated orbit which is not attracting (from the 
numerical simulations it seems to happen for K > 5.11). For larger K the 
values of c n may not decay as fast so we may need to diagonalize the first 
coefficients of Fi as mentioned in the proof. 

Of course sin is periodic so there exist infintely many such orbit that 
differs by (k G Z). As mentioned for K > 5.11 is seems that these 
periodic orbits stops to be attracting and orbits from numerical simulations 
jump ±27r from one periodic orbit to another. This suggests heteroclinic 
connections and chaos but proving it would require some new ideas. 
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